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Abstract 

In the preceding papers (see [1, 2]), the superfluidity of the classical liquid was 
proved under the assumption that the parameters N and r, where N is the particle 
number and r it the capillary radius, tend respectively to infinity and to zero so 
that jj <C where R is the capillary length. In the present paper, this assumption 
is removed. 



1. We first note that solutions of the variational equation for the Vlasov equation 
do not coincide with the classical limit for the variational equations corresponding to the 
mean-field equations in quantum theory. We consider mean-field equations of the form 



ih%-<p\x) = ( -±L-A + W t (x))<p\x), W t (x) = U(x) + I V{x,y)\^{y)\ 2 dy, (1) 



dt \ 2m 

with the initial condition <p\ t =o = (po, where (p G W / 2 00 (R I/ ), J dx\ip (x)\ 2 = 1. 

To find an asymptotic representation of the complex-germ type [3], we consider a 
system consisting of the Hartree equation and its conjugate equation. We then take the 
system of variational equations for it and replace the variations 5tp and 5ip* with the 
independent functions F and G. For F and G, we obtain the system of equations: 

dFHx) t ( 5 2 H . 5 2 H A 

l — = J dy WixWyf iy) + wvfe) 6 iy) J ; (2) 

dGHx) f , ( 5 2 H „ f/ , 5 2 H 



Of 



( S 2 H pt + 6 2 H Qt A 
\6<p(x)5<p(y) 5ip(x)8*ip(y) J 



Roughly speaking, the classical equations can be obtained from the quantum ones using 
a substitution of the form ip = x e ^ S (the WKB method), ip* = x* e ^ S * \ S — S*, x — 
X(x,t) eC°°, S = S{x,t) E 

For variational equations, it is natural to vary not only the limit equation for x an d 
X* but also the functions S and S*. This gives an important new term in the solution 
of the equation for collective oscillations. We consider this fact for the simplest example 
investigated in the famous work by Bogoliubov on a weakly nonideal Bose gas [4]. 

Let U = for (1) in a three-dimensional box with edge length L, and let the L- 
periodicity condition be imposed on the wave functions in this case (i.e., a problem on 
the torus with generators L, L, and L is considered). Then the function 

<p(x) = L- 3/2 e l/h{px - nt \ (3) 
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with p = 2irn/L, where n is an integer-valued vector, satisfies (1) with 

n = ^ + L~ 3 [ dxVix). (4) 
2m J 

We consider functions F^ x \x) and G^(x), where A = 2irn/L,n ^ of the form 

pW\ x ) = L- 3/2 p x e^ p+x)x+il3 ~~ n)tl , 

GW'(i) = L - 3 / 2 a x e^- p+x ^ +n ^, (5) 



where 



H 2 - W 2 = l, Vi = L- 3 [ dxV(x)ei x *. 



From system (6) we obtain 




& = -pA + W — + -V A 2 . (7) 



In this example, we have u = e^ s ^ and w* = e ( h' ) , where s(x, t, ) = px + and 
the variation of the action for the vector (Su, 5u*) is equal to \x ± Qt. 
A more thorough passage to the limit gives V\ — > Vo = L -3 / dxy(x). 
In the classical limit, we thus obtain the famous Bogoliubov relation (7). In this case, 
we have u(x) = 0, and the exact solution coincides with the classical one as in the linear 
Schrodinger equation. The situation with u(x) ^ was investigated in [5], and it turns 
out that a relation similar to (7) is the classical limit as h — > for the variational equation 
in this general case. The curve for the dependence of (3\ on A is called the Landau curve, 
and it specifies the superfluid state. The value A cr at which the superfluidity disappears is 
called the Landau criterion. Bogoliubov explains the superfluidity phenomenon as follows: 
"the 'degenerate condensate' can move without friction relative to elementary excitations 
with an arbitrary sufficiently small velocity" ([4], p. 210). 

But this mathematical consideration is not related to the Bose-Einstein condensate; 
merely the quasi-particle spectrum determined for A < A cr is positive. This means that it 
is metastable (see [6]). The Bose-Einstein condensate is mentioned here only to disprove 
the idea that it follows from what was said above that this consideration applies to a 
classical liquid. 

Indeed, for example, the molecules of a classical undischarged liquid are Bose particles 
if the number of neutrons in the molecule is even. Because every particle (molecule) is 
neutral and is formed of an even number / of neutrons, an iV-particle equation can be 
written for this liquid. Thus, every its particle is 3(2k + Z)-dimensional, where k is the 
number of electrons; there is a dependence on the potential u(xj), Xi G _R 6fe + 3 '; and an 
equation for iV particles x iy i = 1, . . . , N with a pair interaction potential V(xi — Xj) can 
be considered. 

But Bogoliubov found only one series for the spectrum of the many-particle problem. 
As Landau wrote, "N. N. Bogoliubov recently managed to find the general form of the 
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energy spectrum for the Bose-Einstein gas with a weak interaction between the particles 
using a clever application of second quantization" ([7], p. 43). But this series is not unique, 
i.e., it does not exhaust the whole energy spectrum. 

In 2001, we suggested the ultrasecond quantization method [8] (also see [9, 10, 11, 
12, 13]). The ultrasecond-quantized Schrodinger equations, like the second-quantized 
ones, represent the iV-particle Schrodinger equation, and this means that the ultrasecond- 
quantized equation is essentially identical to the original iV-particle equation: it coincides 
with the latter on the 3iV-dimensional space. But in contrast to the second-quantized 
case, replacing the creation and annihilation operators with c-numbers does not yet give 
the correct asymptotic representation; it turns out that its results coincide with those 
obtained by applying the Schroder variational principle or the Bogoliubov variational 
method. 

For the Bardeen exotic potential, the correct asymptotic solution coincides with the 
one resulting from applying the abovementioned ultrasecond quantization method. For 
potentials of general form, in the case of pair interaction for example, the answer turns out 
different. In particular, the ultrasecond quantization method gives some other asymptotic 
series of eigenvalues corresponding to the iV-particle Schrodinger equation, which, in 
contrast to Bogoliubov series (7), are not metastable. They correspond to vortex filaments 
[14]. 

It turns out that the decisive factor here is not the Bose-Einstein condensate but the 
thickness of the capillary (nanotube) in which the liquid flows. If we consider a liquid in 
a capillary or nanotube, then the velocity corresponding to metastable states is not small 
for a sufficiently small radius. Consequently, the liquid flows without friction for a smaller 
velocity. 

The no-flow condition on the boundary of the nanotube (absence of flow) is the Dirich- 
let boundary condition or the Born-von Karman boundary condition. It generates a 
standing wave that can be interpreted as a particle-antiparticle pair: a particle with the 
momentum p orthogonal to the tube wall and an antiparticle with the momentum —p. 

In the boson case, we consider a short-range interaction potential V(xi — Xj). This 
means that only interaction with finitely many particles is possible as N — > oo (N is the 
number of particles). Consequently, the potential depends on N as V N — V^Xi—x^N 1 ^). 
If V(y) is finitely supported in Qy, then the number of particles captured by the support 
is independent of N as N — > oo. As result, superfluidity occurs for velocities less than 
min(A cr , where R is the nanotube radius. The upper bound is determined by the 

condition that the radius of action of the molecule must be less than the radius of the 
nanotube. 

We now present our own considerations that do not relate to the mathematical presen- 
tation. Viscosity is connected with collisions of particles: the higher the temperature is, 
the greater the number of collisions. In a nanotube, there are few collisions because only 
those with the tube walls occur, which is taken into account by the series obtained below. 
Precisely this fact rather than the presence of the Bose condensate leads to the weakening 
of viscosity and consequently to superfluidity. In other words, even for liquid He4, the 
main factor in the superfluidity phenomenon is not the condensate but the presence of a 
thin capillary [15, 16]. 

2. In this part of the paper we refine solutions of equations in variations (26)-(28) and 
(35)-(37) presented in the paper [18]; see also [17] and [19]. 

For the Fermi liquid (for instance, for helium 3), we solve the mathematical problem 
on the reduction of the iV-particle Schrodinger equation with the pair interaction poten- 
tial typical for helium: repulsion as a pair of particles approaches and attraction as the 
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pair of particles moves away. This means that we do not take into account the possible 
resonance interaction of Cooper pairs, as is usual in supercondactivity problems. More- 
over, our assumption forbids the radius of the capillary to be less than the radius of the 
molecule itself (otherwise the molecule simply cannot enter the capillary). However, as 
the radius decreases, a superfluidity domain can occur for a liquid with an even number of 
neutrons. Nevertheless, for a Fermi liquid, superfluidity also occurs, but in quasithermo- 
dynamics rather than in the thermodynamical limit. The notion of quasithermodynamics 
was recently introduced by the author [20]. 

3. The eigenvalues of the equations in variations for Bose-liquid are [20]: 



Ai 



h 



kik 2 ,i = h(k 2 + + \ 



ik 2 ,l + yj H 2 ,i ~ ^Vk 2 ,l 



(8) 



where 



6w = « 2 ((*? - klf + (l 2 - kl) 2 ) + 

+ a {l\ (vi +3k2 - v 2k2 ) + I 2 (uj-fca - v 2k2 ) - k 2 2 (vi_ k2 + v t+3k2 - 2v 2k2 )) - 

- (v l+k2 + v 2k2 )(v l+3k2 + Vi_ k2 - 2v 2k2 )/2, 

Vk2 j = a(2a {k\ - k 2 2 {l\ + I 2 ) + III 2 ) - {l\ + l 2 - 2k 2 2 ) (v l+k2 + v^- 
■ (2a 2 (k 4 2 - k 2 2 {l\ + I 2 ) + l 2 J 2 ) + 

+ a{l\ (2v t - k2 + v i+k2 - v 2k2 ) + I 2 (2v l+3k2 + v l+k2 - v 2k2 ) - 

- 2k\ {vi +3k2 + vi_ k2 + v l+k2 - v 2k2 ) ) + 

+ 2(v i+3k2 - v 2k2 )(vi- k2 - v 2k2 ) + (v l+k2 + v 2k2 )(v l+3k2 + vi-fr - 2v 2k2 ) s j/A] 
and where the notation 

h 2 

a = - — , h — I + 2k 2 . 
2m 

was used for brevity. 

Relation (8), as compared with (27) in [17], is uniform with respect to k 2 as k$ 
Note that, if k 2 = 0, then the Bogolyubov relation holds, 



M,k lt i — k\l + 

m 



2m 



+- vi ] - vf. 



oo. 



(9) 



For a system of identical Fermi-particles, the eigenvalue problem for the system of 
equations in variations can similarly be reduced to the problem of finding the eigenvalues 
of the equation 

AX = MX. (10) 
Here A = A + ^k±(k 2 + I), X - is a column vector of the form 



X 



u 2) i 



M - stands for the matrix 
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with the elements 



B± = B k2 j + 



v i-k 2 



2 ' 



B 2 — B k2j i + 2k 2 + 
Mi = 2i(u,_ fca - u )<Pk2,j, 



Vl 

V 2 



Vl+k 2 ~ v 2k 2 



Vl-k 2 ~ v l+k 2 



v l+3k 2 ~ v l+k 2 



(11) 



M 2 = 2i(v - v l+3k2 )Lp k2t i +2k2 , F = i(v 2k2 - v i+k2 ) (<£ k2) i +2k2 - ip k2t i) , 
where the numbers B k2 j and ip k2 j are 



£ 2 



2 6, 



*2>1 , ^ h 2 



2 ± yrw _1 ' 



(I 2 - kl) + (v - v 2k2 ) 



(12) 



^_ fc2 - Vi +k2 



fci, k 2 and Z are three-dimensional vectors of the form 

2n (— — — ^ 

VLi L2 L2/ 



(13) 



ni, n2 and n 3 - n3 are integers. The summation is taken over all values of ni, n 2 , n 3 , and 
lk 2 , = c?({ll-kl) 2 + {l 2 -kl) 2 ) + 

+ a {l\ (v l+3k , 2 - v 2k2 ) + I 2 {vi_ k2 - v 2k2 ) - k\ (vi_ k2 + v l+3k , 2 - 2v 2k2 )) + 
+ (vi+k 2 - v 2 k 2 )(vi-k 2 + vi+3k 2 - 2v 2k2 )/2, 

tj k2 ,i = a(2a {k\ - k 2 2 {l\ + I 2 ) + l\l 2 ) + {l\ + l 2 - 2k 2 2 ) (v l+k2 - v 2k2 ) ) • 
• (2a 2 (kt - k 2 2 (l\ + I 2 ) + l 2 ,l 2 ) + 

+ a(ll (2vi- k2 - v 2k2 - vi +k2 ) + I 2 (2vi +3k2 - v 2k2 - vi +k2 ) - 
- 2kl (vi +3k2 + vi- k2 - vi +k2 - v 2k2 )) + 

+ 2(v i+3k2 - v 2k2 )(vi-k 2 - v 2k2 ) - (v l+k2 - v 2k2 )(vi+3 k2 + Vi- k2 - 2v 2fc2 ))/4, 
where the notation 

a = - — , li — I + 2k 2 . 
2m 

is used for brevity. The eigenvalues of the system of equations in variations are 



fv 



Ai.kifca.i = h(k 2 + + \ 

m \ 



Zk2,l + \/ 'ffc 2 ,I - ^Vk 2 ,l 



(14) 



For fc 2 = 0, i.e., along a capillary, the superfluidity in the quasithermodynamics holds and, 
for &2 7^ 0, i.e., if there is a reflection from the walls of the capillary, then an instability 
occurs because the eigenvalues become complex by (14). This gives the same upper bound 
for the critical speed as that in the Bose-case (due to occurrence of a vortex). 

According to [17, 18, 19] only the thermodynamical and quasithermo-dynamical limits 
in statistical physics exist. The above solution on the superfluidity of a Fermi liquid 
enables us to claim that the phase transition in this problem from the superfluid state to 
the normal one is an example of a phase transition in quasithermodynamics. 

The results obtained in this paper are based on the author's paper of the year 1995 
[5] which is given in appendix below. 
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Quasi-Particles Associated with Lagrangian Manifolds 
Corresponding to Semiclassical Self-Consistent Fields. 
Ill 1 

In the preceding part of this paper, we presented Eqs. (25) for quasi-particles associated 
with an n-dimensional Lagrangian manifold and Eq. (29) for quasi-particles corresponding 
to a (2n — l)-dimensional manifold. These equations were written out only in the x-chart, 
and the quantum corrections were given without proof. In this part we essentially use 
the canonical operator method to obtain Eq. (25) with corrections in the x-chart as well 
as in any other chart of the canonical atlas [1]. To derive the correction in Eq. (29), a 
"modified" 5-function must be used, and this will be done in the next part of the paper. 

To obtain the result in an arbitrary canonical chart, one should pass on to the p- 
representation with respect to some of the coordinates in the Hartree equation. This is 
actually equivalent [2] to considering the Hartree-type equation 



i 



H (x,-ih-^^) + J dyi;*(y)H 1 (^x,-ih-^;y,-ih^^(y)^(x) = Q^(x), (A.l) 

where x,y G M n , ip G L 2 (R n ) is a complex-valued function, h > 0, Q G R, and the indices 
1 and 2 specify the ordering of the operators x and —ihd/dx. The function Hi satisfies 
the condition Hi(x,p x ;y,p y ) = Hi(y,p y ;x,p x ). Equation (A.l) generalizes the ordinary 
Hartree equation (Eq.(l) in [4], where N = 1). The study of Eq. (A.l) is important, for 
example, if one makes an attempt to find a solution to the Hartree equation (1) in the 
momentum representation, 



(i/h)px 



dp 



(27r£)™/ 2 ' 



Let us also discuss the variational system associated with Eq. (A.l), which can be 
obtained as follows. Along with Eq. (A.l), let us write out the conjugate equation and 
consider the variations of both equations assuming that the variations 5ip = F and 5ip* = 
G are independent. 



^.P.Maslov, Russian J. Math. Phys. 1995, v.3, N3, 401-406. 



The variational system has the form 



F(x) 



1 11 
H (x,-ih—^J -n + J dy^^H^x,- ih—;y,- ih—^ij(y) 

i i 
+ J dyfG^H^x,- ih—\y,- ih—^ip(y) 
i i 

+ r{y)Hi ( 1 - ih^-- y, - ih^- )F( 2 /))^(x) = -PF(x), (A.2) 

2 2 2 

H (x,ih— / dyij^H^xJh—^yJh—^ij^y) G(x) 



+ 



dy 

2 

dy[F{y)H 1 ( K h,ih—-y,ih— ^*(y) 

2 2 



Equations (A.l) and (A.2) play an important role in the problem of constructing asymp- 
totic solutions to the A^-particle Schrodinger equation as N — > oo [5]- [7]. 

For example, the spectrum of system (A.2) (possible values of (3) corresponds to the 
spectrum of quasi-particles. Namely, the difference between the energy of an excited 
state and the ground state energy is given by the expression J2 k 0k n k, where the numbers 
rife G Z + , k — 1, oo, which are equal to zero starting from some k, define the eigenfunction 
and the eigenvalue of the excited state, and /3k G K are the eigenvalues of system (A.2). 

In this paper we are interested in asymptotic solutions to Eqs. (A.l) and (A.2) as the 
"inner" h tends to zero. 

Asymptotic solutions to Eq. (A.l) are given [8] by the canonical operator on a La- 
grangian manifold A n = {x = X(a),p = P(a)} invariant with respect to the Hamiltonian 
system 

dH(x.p) • dH(x,p) .. „ N 



where 



H(x,p) = H (x,p) + J dn a Hi(x,p;X(a),P(a)), 



a G A™, and dfi a is an invariant measure on A™. The Lagrangian manifold lies on the 
surface H(x,p) = Q. If a chart A is projected diffeomorphically in the x-plane, then the 
canonical operator acts as the multiplication by exp{(i/h)S(x)} /y/1, where S(x) = J pdx 
on A n and J = Dx/Dfi a . We are interested in finding asymptotic solutions to Eqs. (A.2). 
Without loss of generality, we can confine ourselves to the case of rc-chart. Indeed, to 
obtain similar expressions in the p-chart, one must consider the Fourier transformation 
of Eqs. (A.l) and (A.2) and apply the same technique, since the form of the equations 
remains unchanged. 

Let us seek the asymptotic solutions to Eqs. (A.2) in the x-chart in the form 

F(x) = f(x)i>(x), G(x) = g(x)r(x), (A.4) 

where the functions / and g, in contrast to tp and ip*, have a limit as h — > 0. One can 
consider a more general case, by allowing / and g to be functions of x and —ihd/dx, but 



in the leading term as h — > we have 

and so we arrive at functions / and g that depend only on x. 

The second equation in system (A. 2) can be rewritten in the form 

H °( X,ih lL) + J dy ^^ Hl { x ' ih ^ y ' ih ^)^*^ y y^(x)^*(x) 



+ 



J dy\i){y)c{y)H l ( K x,ih^y,ih^^*(y) (A.5) 

Is 

dy, 



+ ^(y)[Hi(x,ih-^-;y,ih-^-y,g(y) i/j* (y)}ip* (x) = (3g(x)ijj*(x), 



where [A; B] = AB — BA and 

c(x) = f(x)+g(x). (A.6) 

Equation (A.l) is used in the derivation of Eq. (A.5). We observe that all terms containing 
the function g on the left-hand side in Eq. (A.l) are 0(h), since the commutator of two 
operators depending on x and —ihd/dx is equal, in the classical limit, to (—ih) times the 
Poisson bracket of the corresponding classical quantities. 

Thus, the function c, as well as the eigenvalue f3, is assumed to be 0(h). Let us rescale 
these quantities as follows: 

c(x) = hc(x), (3 = hp. (A. 7) 

Now we can derive the equation for g,c, and (5 in the leading term in h and the first 
correction to it from Eq. (A.5), making use of the following relations: 

i) \ A ( x ,ih— = Yih—(k,ih—)^--Y h " 82 A (kih- ) d ^ 

L V dx J ' J ^—f dp a V dx J dx a ^ 2 dpadpb V dx / dx a dxb ' 

a=l ' a, 6=1 

(A.8) 

where p a = ihd/dx a , A(x,p) is a function R 2n — > C, £ : IR n — > R; 

ii) V'Oe) — x( x ) h)e^/ h ^ s ^ x \ where x — l/V~J in the leading term in h; 

Hi) th^-mm*) = e -m)s( X) ( 9s + ^ 

ax V or or / 



dx dx J \dx/ , dp a dx a 2 ^ dp a dpbdx a dxb 

a=l a, 6=1 

, (^) 2 y d*B g | (^) 2 ^ a 2 g g 

2 ^dp a dp b dx a dx b 2 j-^ dp a dp b dp c dx a dx b dx c 

a, 0=1 a,6,c=l 

(^) 2 ,A <9 3 5 <9 3 ,S 



6 j-^ dp a dp b dp c dx a dx b dx c 

a,b,c=L 



8 fc ^ i dPadpbdp c dp d dx a dx b dx c dx d 

a,b,c,a=l 

where all derivatives of B are evaluated at the point p = dS/dx. 



in 



These relations can easily be obtained for monomial functions A and B. An application 
of formulas i)-iv) yields the equation 



dpi dx a 

a=l 



.. ,,dH x h dg , w .. d 2 H <91nJ\ . 



dpi 1 
d 3 H 



a,b=l 

d 2 S 



dp x a dp x dx b 



a,b,c=l 



dp%dpfdp x dx b dx ( 



<*<<»» - I E i^rPW) 



+ 



Y 

~2~ 



n 



Z_-^ Vdp x a dpl dx a dx b 



(X(a)) 



2 ^ dx a dx b 

a, b=l 

d 2 H l d 2 S 
dpldpl dx a dx b 



dp x a dp. 



, r9ffi91nJ.^. ,, dHid\nJ.. x , . 

d^c(X(/3)) E | (*(«)) + ^5T — (^(^) 



+ 



t/^E|fw»{E 



d 2 H x d\nJ, v ,„„ d 2 H l din J 



d 3 H t 



d 2 S 



y , 

6 '~ L ^ dpldpl dp > x dx b dx, 



hf d y d 2 H x d 2 g 
2 J P dpldp y b dx a dx b 



(X(a)) + 



d 3 H 1 d 2 S 



dpldpl dx b v " vr "" 1 dpldpl dx b 



(*(«))) 



dpldp y b dp v c dy b dy c 



(X((3)))} 



(X(P)) = 0; 



in this formula the arguments 

x = A», p* = P(cO, y = X(/3), p" = P(/3) 



(A.10) 



(A.11) 



of the function Hi and of its derivatives, as well as the arguments x = X(a), p x = P{a) 
of the function H, are omitted. 

Let us now find another equation relating g to c. To this end, let us multiply the first 
equation in system (A. 2) by ip*(x) and the second equation by ip(x). Let us subtract the 
first product from the second. We obtain 

d 



(3ip*(x)if;(x)c(x) = ip(x) H(x,ih—y,g(x) ip*(x) 

+ ip*(x) H^x,—ih—^;g(x) ip(x) — ip*(x) H^x, —ih—^;c(x) ip(x 
+ i){x) J dyip(y) 



i d d \ i 



i/j*(x) I dy^*(y) g(y);Hi(x : -ih^-;y : -ih^j ip(y)4>(x) 



d 



(A.12) 



+ i/>(x) J dyip(y)c{y)Hi^x,ih^;y,ih 



d_ 

dy 



dy< 



-ip*(x) J dy7p\y)Hi(x,-ih^;y,-ih^jc(y)'ip(y)'il;(x). 

Let us use Eqs. (A.8)-(A.10). We find the following equation for g and c modulo 
0{h 2 ): 



d 



1 1 



*tZ^( X < a »- - t £L SI (Ma)) 



dp x a dx a 
^ d 2 H 



^^P? dx a dx b 



a,b= 



\ dPadpl dx a 



dx b 



_ y f X (a)) 9 " H 

^ dx a dpldpldp? dx b dx c 

d 2 S 



a,b,c=l 



d 2 9 



d 2 E x 



^ + ^~ b {Xm dpldpldpl) 



a,b=l 

d 2 S 



d 3 H l 



+ 



d 2 H l d In J 



1 / ^c(X(/3))JJ(^— (X(a))^ + ^— (X(/3)) 
0=1 



dpldp y b 



dx a 



<9pS <9a; a 



a=l 



+ i 



^ 02 77 o2 o 

^c(X(/5)) £ f ^JL _ _ ( X (a)) 



a,6,c=l 



<9c 



^ \dp x a dp x b dx a dx b 
d 3 H d 2 S 



d 2 H x d 2 S 
dpadpl dx a dxb 



W))) 



dPadpldp* dx b dx ( 



(X(a)) 



h ^ «9c <9 2 iJ 9 In J /i >A <9 2 iJ <9 2 c 

{a)) dpidpi dx b ( + 2 ^ dpidp x b dx a dx b [ {a)) 



a,b=l 



h f dc d 3 H 1 

ft ,0 ; C — -L 

+ 2j d ^^[dpJdpj{dx~a 



d 2 S /v . d 3 ^ 



<9 2 S 



ids*. 



*(X(0)^(X(a))** 



""" (.Y(;3)) - ^-(X(ffl) 



dpldpfdpc dx b dx, 
d In J 



8x n 



dx h 



dx a 



dx h 



dpldp] 



0. 



(A.13) 



If H (x,p x ) = pl/2 + U(x) and Hi(x,p x ;y,p y ) = V(x,y), then Eqs. (A. 10) and (A.13) 
become much simpler and acquire the form 



(iVSV-(3)g + J V(x,X(a'))c(X(a'))dii a! + ^(-Ag + V In JVg) = 0, 
(iVSV - /?)c - Ag + V In JVg - ^(-Ac + V In JVc) = 0. 



(A.14) 



From Eqs. (A.14) one can approximately find the functions F and G, which are important 
for constructing approximate wave functions in the iV-particle problem as iV — > 00 [5]. 

Let us now relate the obtained results to the solution to variational equation for the 
Vlasov equation, obtained in the preceding part of this paper [3]. 
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Let p be the projection on the function ip. Its kernel is p(x,y) = ip(x)ip*(y), and its 
symbol is p(x,p) = ip(x)^p* (p)e^^^ px . The operator p satisfies the Wigner equation, which 
reduces to the Vlasov equation as h — > 0. The operator a with the kernel F(x)ip*(y) + 
ip(x)G(y) is equal to 

S = fp + P9 (A.15) 

and satisfies the variational equation to the Wigner equation, which is reduced to the 
variational equation for the Vlasov equation (20) obtained in [3]. In Eq. (A.15) / and 
g are the operators of multiplication by the functions / and g. We see that in the 
semiclassical approximation the symbol of a is 0(h), since a = [p;g] + he p and 

a(x, P )c,h(-i£^(x, P )^ + cp). 

Since p is the (^-function in the semiclassical approximation [3], the function a is actu- 
ally the sum of the (^-function and its derivative. Equations (A. 14) are consistent with 
Eqs. (24) obtained in [3] for the coefficients of 5 and 5'. Thus, the approach suggested in 
this part allows us to find an asymptotic formula for er as well. 

The author is deeply grateful to O. Yu. Shvedov, whose assistance in carrying out all 
computations was invaluable. 
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